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■ ABSTRACT 



Aims. We have studied the molecular content of the circumstellar environs of symbiotic stellar systems, in particular of the well know 
objects R Aqr and CH Cyg. The study of molecules in these stars will help to understand the properties of the very inner shells around 
the cool stellar component, from which molecular emission is expected to come. 

Methods. We have performed mm-wave observations with the IRAM 30m telescope of the I2 CO 7=1-0 and 7=2-1, l3 CO 7=1-0 
and 7=2-1, and SiO 7=5-4 transitions in the symbiotic stars R Aqr, CH Cyg, and HM Sge. The data were analyzed by means of 
a simple analytical description of the general properties of molecular emission from the inner shells around the cool star. Numerical 
calculations of the expected line profiles, taking into account the level population and radiative transfer under such conditions, were 
also performed. 
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Results. Weak emission of 12 CO 7=1-0 and 7=2-1 was detected in R Aqr and CH Cyg; a good line profile of 12 CO 7=2-1 
O . R Aqr was obtained. The intensities and profile shapes of the detected lines are compatible with emission coming from a very 

small shell around the Mira-type star, with a radius comparable to or slightly smaller than the distance to the hot dwarf companion, 
10 14 - 2 x 10' 4 cm. We argue that other possible explanations are improbable. This region probably shows properties similar to 
.£3 \ those characteristic of the inner shells around standard AGB stars: outwards expansion at about 5-25 km S" 1 , with a significant 

acceleration of the gas, temperatures decreasing with radius between about 1000 and 500 K, and densities ~ 10 9 - 3 x 10 8 cm 4 . Our 
model calculations are able to explain the asymmetric line shape observed in I2 CO 7=2-1 from R Aqr, in which the relatively weaker 
red part of the profile would result from selfabsorption by the outer layers (in the presence of a velocity increase and a temperature 
decrease with radius). The mass-loss rates are somewhat larger than in standard AGB stars, as often happens for symbiotic systems. 
In R Aqr, we find that the total mass of the CO emitting region is ~ 2 - 3 x 10~ 5 M , corresponding to M ~ 5 x 10~ 6 - 10~ 5 M Q yr~' , 

■ and compatible with results obtained from dust emission. Taking into account other existing data on molecular emission, we suggest 
' that the small extent of the molecule-rich gas in symbiotic systems is mainly due to molecule photodissociation by the radiation of 

the hot dwarf star. 

Key words, radio lines: stars - stars: circumstellar matter, mass-loss - stars: binaries: symbiotic - stars: individual: R Aqr, CH Cyg 
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■ 1. Introduction particularly important. However, molecules are abundant only 

in the innermost regions of the CSEs around SSs. Bands of 

Symbiotic stellar systems (SSs) are very close binary systems CQ TiQ and other molecu les are often detected in SSs, com- 
. . composed by a cold red giant and a very hot dwarf compan- ing from ±e outer atmospher e f the cool primary, as in stan- 
JX; ■ ion. The strong interaction between both stars yields a num- dard red giants However, to our knowledge, molecular low- and 
H i ber of interesting and sometimes spectacular phenomena, like intermediate-excitation emission (SiO masers, CO thermal lines, 
«i , the ejection of fast colhmated flows and the formation of high- etc)> known to come from the circumste llar shells, ha d been only 
excitation, bipolar nebulae (see e.g.| Corradi et al.| | 2003[ and ref- detected in two SSs: R Aqr and Hl-36 (llvison et al]ll994 ll998. 
erences therein). They are thus a very attractive laboratory for ISchwarz et all [1991 ISeaauist et al.lll995l) . SiO maser emission 



studying various aspects of stellar evolution in binary systems (y>0 lines at mm wave i engths ) in R Aqr is relativel y 'normal', 

and of the circumstellar chemistry and structure under these (ex- com pared to that o bserve d in other AGB stars (e.g. iPardo et al.1 

treme) conditions. [2004 ICotton et al.ll200H iKamohara etal]|2010l) . We recall that 

Some symbiotic systems show strong IR excesses, due to the SiO maser emission comes from the innermost circumstellar 

emission of dust grains formed in the circumstellar envelopes reg ions, in our case at about 5 x 10 13 cm. H 2 masers have 

(CSEs) ejected by the cold primary. In these objects, the ejec- been also detected in these same two objects, with character- 

tion of circumstellar gas from the primary is thought to be i SI ics that suggest that H 2 Q emission forms only a t a few stellar 

radii ([Seaquist et al. 1995, I vison et al Jl994ll998h .InRAar.an 

* Based on observations carried out with the IRAM Pico Veleta H^O-rich shell extending also ~ 5 X 10 13 cm has been detected 

30m telescope. IRAM is supported by INSU/CNRS (France), MPG from observations of H 2 vibrational bands (Ra riand et~aT] 
(Germany) and IGN (Spain). 
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2008) . Molecules characteristic of the outer shells are yet weaker 
and very rarel y detected. OH masers have been wel l detected 
only in HI -36 dlvison et al.lll994ISeaquist et al.lll995l) . The CO 
thermal lines, which in general come from still outer shells (of- 
ten within about 10 17 cm), are extremely weak in SSs. Previous 
to our wor k, only a very tentative det ection of R Aqr had been 
reported bv lGroenewegen et al.1 ([199 91). In fact, our observations 
(see next sections) show that the feature detected in R Aqr by 
these authors is mostly due to baseline ripples in the spectrum, 
since the actual intensity of the line is about three to four times 
lower. 

It seems well established that the detectability of molecu- 
lar emission in SSs increases for lines coming from very com- 
pact regions around the AGB star and for s ystems showing rel- 
atively large dista nces between the stars dSchwarz et al.l 1995, 
llvison et al1ll998l) . Therefore, the lack of detections of molecu- 
lar lines is probably due to photodissociation by the UV radia- 
tion from the hot companion or to dynamical disruption of the 
emitting regions (see further discussion in Sect. 5). 

The properties of the innermost shells around red giants, 
within about 10 stellar radii, from which molecular lines in SSs 
seem to come, are not very well known even for i solated stars. 
Both theoretical and observational studi es (e. g . IHinkle et al.l 
[1981 lHofner etal.1 [l998l lAndersen et all I2IM ISandinf feoOS) 
suggest that relatively high temperatures, between 500 and 1000 
K, and densities, 10 8 - 10 9 cm -3 , are present. These inner shells 
are thought not to show the fast expansion characteristic of the 
outer regions. They are probably pulsating, due to shocks origi- 
nated from the photospheric pulses, or show incipient expansion, 
since in these regions dust grains are being formed and radiation 
pressure efficiently acts onto them. In symbiotic systems, pulsa- 
tion and outwards acceleration may also dominate the dynamics 
of the inner circumstellar layers, but the gravitational effects of 
the secondary cannot be neglected, since SSs are interacting sys- 
tems. In particular, it is remarkable that SSs present mass-loss 
rates systematically larger than those of isolated AGB stars (e.g. 
Mikolaiewska 1999, and references therein), which may be due 
to such gravitational effects. 

In this paper, we present observations of molecular mm- 
wave lines in three SSs: R Aqr, CH Cyg, and HM Sge. 12 CO 
emission is well detected in R Aqr and CH Cyg. 

R Aqr and CH Cyg are bright and nearby symbiotic stars, 
extensively studied over the whole spectral range. The dis- 
tance to R Aqr has been accurately determined from recent 
VLBI measurements of its parallax by Kamohara et al. (2010); 
we will adopt their distance value, D = 214+ 4 2 i . A distance 
of 244 +49 f pc was m easured for CH Cyg from Hipparcos data 
(Ivan Leeuwenl l2007h . In both systems, the hot component is 
an accreting white dwarf showing spectacular activity: irregu- 
lar accretion-powered outbursts accompanied by massive out- 
flows and jets (e.g. Kellogg et al. 2007, Karovska et al. 2007, 
and references therein). The cool component in R Aqr is a Mira 
variable with a pulsation period of 387 d , whereas in CH Cyg 
it is an M7III semiregular variable with a complex variabil- 
ity (e.g. Gromadzki & Mikolaiewska 2009, Mikolajewski et al. 
1992, and references the rein). Both have relatively well 
known orbital par ameters ((Gromadzki & Mikolaiewska 2009, 
IHinkle et alll2009l) . Aqr and CH Cyg have the longest orbital 
periods measured in well studied symbiotic systems, respec- 
tively 43.6 yr and 15.6 yr. The o rbital solution for R Aqr 
(iGromadzki & Mikolaiewska 120091) implies that the average 
component separation is ~ 2.25 x 10 14 cm (15 AU). During our 
observations the component separation was ~ 17.6 and 16.5 AU, 
in May 2008 and May 2009, respectively. In the case of CH Cyg, 



the orbital elements from Hinkle et al. (2009) yield an average 
component separation of ~ 9 AU (1.35 x 10 14 cm), and during 
the May 2009 observation the separation was ~ 9.8 AU. The 
red giant radius in both systems is known from interferometric 
measurements, ~ 1 .9 AU in R Aqr (Gromadzki & Mikolajewska 
2009, and references therein), and ~ 1.2 AU in CH Cyg (e.g. 
Dyck et al. 1998), with uncertainty set mostly by the uncertainty 
in their distances, ~ 15 % in both cases. 

The symbiotic nova HM Sge is composed of a Mira variable 
with a pulsation period of 527 d , embedded in an optically thick 
dust shell, and a white dwarf companion slowly declining from 
a thermonuclear nova outburst started in 1975 (Belczynski et al. 
2000, and references therein; Muerset & Nussbaumer 1994). The 
orbital period is unknown, likely higher than ~ 100 yr. Eyres et 
al. (2001) measured the binary component positions using HST 
images, and estimated a projected binary separation of 40 ± 9 
mas, and a position angle of the binary axis of 130 + 10 de- 
grees, in agreement with that suggested by Schmid et al. (2000) 
based on spectropolarimetry. Unfortunately, the distance to HM 
Sge is rather uncertain - published values range from 0.3 to 3.2 
kpc (e.g. Richards et al. 1999). Eyres et al. adopted D = 1.25 
kpc, and a component separation of 50 AU. The recently re- 
vised period-luminosity relation for single Miras (Whitelock et 
al. 2008) would place HM Sge at D = 2.5 kpc; this estimate is 
however uncertain due to the peculiar nature of the object, par- 
ticularly because of the low amplitude and poor periodicity of its 
light curve (see the AAVSO database). 

As we will see, the CO lines in these systems are extremely 
weak, typically about 100 times weaker than for standard AGB 
stars. We will argue that this low intensity, as well as the other 
main properties of the detected lines, show that the CO emis- 
sion only comes from very inner circumstellar regions around 
the cool stellar component, closer than about 1-2 10 14 cm. 



2. Observations and data reduction 

We have used the IRAM 30m telescope, at Pico de Veleta 
(Spain), to observe the mm-wave emission of 12 CO and 13 CO, 
J = 1-0 and J = 2-1, and of SiO 7=5-4, in the symbiotic stars 
R Aqr, CH Cyg, and HM Sge. As calibration standards, we also 
observed the sources IRC +10216, CRL2688, and NGC7027, 
whose CO emission is strong and well characterized. 

The observations were performed in two observing sessions, 
in May 2008 and May 2009. During the first session we used 
the A 100 and B100 receivers for the 3 mm band, and the A230 
and B230 for the 1 mm band, to observe four lines simultane- 
ously. The data were recorded using the 1 MHz filterbank and 
the VESPA autocorrelator. For the second session, we used the 
new EMIBQ receiver, observing simultaneously in the E090 and 
E230 bands (3 mm and 1 mm) in dual polarization mode. In this 
run, only the 12 CO lines were observed. The data was recorded 
using the WILMA and VESPA autocorrelators. 

The spatial resolution of the observations at 3 mm is 23", 
and 12" at 1 mm wavelength. Frequent pointing measurements 
were performed to measure and correct pointing errors; typi- 
cally, errors no higher than 3" were found, which practically 
have no effects on the calibration. The observations were done 
by wobbling the subreflector by 2' at a 0.5 Hz rate. This method 
is known to provide very stable and flat spectral baselines 

The atmospheric conditions during the observations were 
good. The average zenith opacity during both observing runs 
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Fig. 1. 12 CO lines detected in R Aqr. 



Fig. 2. I2 CO lines detected in CH Cyg. 



3. Line formation in inner circumstellar shells 
around symbiotic stellar systems 



was - 
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0.2, slightly better at 1 10 GHz and slightly worse at 1 15 31 ■ Sim P le Predictions of molecular line intensity 



The data presented here has been calibrated in units of 
(Rayleigh-Jeans-equivalent) main-beam temperature, corrected 
for the atmospheric attenuation, r m b, using the standard chop- 
per wheel method. Calibration scans (observation of the hot and 
cold loads, and of the blank sky) have been performed typically 
every 15-20 minutes. The temperature scale is set by observing 
hot and cold loads at ambient and liquid nitrogen temperatures. 
Correction for the antenna coupling to the sky and other looses 
have been done using the latest values for these parameters mea- 
sured at the telescope. The accounting for the sky attenuation 
is computed from the values of a weather station, the measure- 
ment of the sky emissivity, and a numerical model for the atmo- 
sphere at Pico de Veleta. Finally, we checked that the calibration 
of the different observations was compatible. We also compared 
the intensities of IRC +10216, CRL2688, and NGC7027 with 
previous observation s to check calibra t ion un certainties; we took 
observational data in lBuiarrabal et al.l ([2001), but taking into ac- 
count that recent measurements suggest that the 1 mm data in 
that paper seem to be overcalibrated by about a 20%. The cor- 
rections applied to the calibration of the different observations 
were always moderate, not larger than ~ 20%, which can be con- 
sidered as a measure of the absolute calibration uncertainty. 

All the data were averaged and rebinned, to get an adequate 
velocity resolution. Baselines of degree 1 or 2 were subtracted. 
We note that in the EMIR observations a source of noise at inter- 
mediate frequencies was detected, in such a way that averaging 
the WILMA and VESPA data a slightly smaller noise was ob- 
tained. (Averages of observations obtained simultaneously with 
different spectrometers is not a standard procedure in radioas- 
tronomy, because the noise from low-frequency parts of the de- 
tection chain is usually negligible and this procedure does not 
help to improve the S/N ratio.) In any case the differences are 
small; as an example, the spectra shown in Fig. 4 (to be com- 
pared to that shown in Fig. 1) was obtained after averaging both 
spectrometers. 



As discussed in Sect. 1, the properties of molecular lines from 
SSs, particularly their weak emission, are thought to be due to 
effects of the companion on the outer CSE, by means of pho- 
todissociation or of dynamical disruption. If the detected CO 
emission comes from regions closer than the separation between 
both stellar components (several AU, 1-2 x 10 14 cm), it must 
present very different properties than the usual circumstellar CO 
lines, which come from regions as far as ~ 10 17 cm. 

The excitation of the low-/ CO transitions is in general very 
easy to describe, the population of the levels only requires tem- 
peratures i 15 K and is thermalized (i.e. in LTE) for densities 
larger than 10 4 cm 4 . Temperatures close to 1000 K and densi- 
ties J 10 8 cm" 3 are expected in the innermost circumstellar lay- 
ers around AGB stars, say at a few stellar radii (Sect. 1). Even 
for extended envelopes, these conditions are sati sfied in practi- 
cally all the emitting region (see discussion in e.g. Tevssi er et al.l 
2006). CO lines from regions at ~ 10 14 cm are expected to be 
very weak compared to those from standard envelopes, due to 
the large dilution factor in that case. 

We can perform simple, but accurate estimates of the emis- 
sion of CO low-/ lines from the inner CSE (the predictions of 
our simple calculations will be checked later on by means of nu- 
merical calculations). Due to the very high densities expected 
in these regions, the populations of all the relevant CO rotational 
levels are accurately thermalized. (It is easy to show that the pop- 
ulation of the vibrationally excited levels is negligible, since the 
vibrational de-excitations are very probable and the stellar IR ra- 
diation is strongly diluted a soon as we are placed at a few stellar 
radii.) Accordingly, the population of each /-level, n(J), is given 
by: 



n(J) = gjx(J)=gj X (0)e- E > /kT t 



(1) 



where x is the population per magnetic sublevel, g is the statisti- 
cal weight (in our case, gj = 2 J + 1), Ej is the energy of level /, 
and TV is the kinetic temperature. Other symbols in our formulae 
(k, h, ...) have their usual meanings. 
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In terms of the total density, n, and the CO relative abun- 
dance, x(coy. 



x(J) = nXe- E ' lkTk /F(T k ) 



(2) 



F(T k ), the partition function, is accurately given in our case (a 
linear thermalized molecule, T k » B) by F(T k ) = T k /B, B being 
the rotational constant in temperature units (for CO, Bh/k = 2.75 
K). We recall that, in our case, for the CO rotational line J — » 
7-1 and properly adapting the units of B: Ej = BJ(J +1) and 
vj - 2B.I. 

With these expressions, we can easily write the representa- 
tive absorption coefficient and optical depth, r, in a medium with 
characteristic kinetic temperature, density, CO abundance, ve- 
locity dispersion (AV) and size (L). We must keep in mind that 
in our case the line profile is given by the velocity dispersion 
within the source. For transition J — > J - 1 : 



t(J) = nX 



1 



F(T k ) Snv* 



C —A jgj [x(J-\)-x(J)]-^ 



AV 



(3) 



The Einstein coefficient in this linear, simple molecule is Aj = 

6 3ft? V^~-> A 4 being the permanent dipole moment that in CO is 
particularly low. Note that here r is the opacity at the central fre- 
quency of the line, v, corresponding to the systemic velocity of 
the source V sys . For nearby frequencies, v + Av/2 (V sys + AV/2), 
r should vary depending on the velocity distribution within the 
source. (Of course, r decreases strongly far from the resonance 
frequency.) We will see below the dependence of r on the ob- 
served velocity (i.e. on the observed Doppler-shifted frequency) 
for the velocity fields expected in our sources, around their sys- 
temic velocities. 

The dependence of t on T k and J is interesting. Since, in 
our case (and in many cases for CO low-/ transitions), we can 
assume that the level populations are thermalized and T k s> B, 



t(J) ~CXn — B A n — 



(4) 



where the constant C only includes usual mathematical and 
physical constants. Accordingly, for different values of the tem- 
perature, 



r(J) 2 1/T k 2 , 

and, for different transitions (and the same 7\), 
t 5 J 2 . 



(5) 



(6) 



The characteristic brightness of a source (neglecting the cos- 
mic background, whose temperature is much smaller than T k in 
our case) is 



B(v) ~ S(l -e- nv >) 



(7) 



Where S , the source function, depends on the frequency much 
less sharply than the line profile: 



2hv i 



x(J) 



c 2 x(J-l)-x(J) 



(8) 



Since we are assuming that the populations are thermalized (eqs. 
1, 2) and, again, that T k » B, we can safely apply the Rayleigh- 
Jeans approximation and write the brightness temperature as: 



T b (v) ~ T k (l - e-™) 



(9) 



Leading to the usual equations for the optically thick case, 

T b ~ T k , (10) 

and for the optically thin case, 

7b(v) ~ T k r(v) , T b (V) ~ T k r(V) . (11) 

Note that equation 10 holds for the range of velocities around 
the systemic one for which the opacity is still higher than 1, and 
that 7b in eq. 1 1 keeps the velocity dependence of t. We also 
note that, if we allow the amount of emitting molecules to vary, 
the optically thick case always represents the upper limit to the 
source emission. 

From these equations we can deduce also the dependence of 
7b on T k , taking in particular into account the dependence of t 
on T k (see eq. 4): 7b °c T k in the optically thick case, and 7b x 
l/T k in the optically thin case. 

Finally, for a source with a typical size L, in arcsecond units, 
observed with a telescope beam with typical half-power beam 
width equal to W b , also in arcseconds, the main beam tempera- 
ture T mb (over the cosmic background) will be: 



T mb ~ T b (L/W b y 



(12) 



[We are assuming that the typical source size L represents the 
source observable diameter at half-maximum, and that the dilu- 
tion factor within the beam, {L/W b ) 2 , is much smaller than one, 
since the source is in our case unresolved.] 

3.2. Line shapes: optically thin emission 

The typical profiles of molecular lines from CSEs around stan- 
dard AGB stars, very extended and expanding at supersonic al- 
most constant velocities, are well known. Parabolic, flat-topped 
or two-horn profiles are expected depending on whether the lines 
are optically thin or thick and on whether the sou rce is spatially 
resolved by the telescope (e.g. lOlofsson et aill982l) . 

The profile shapes are of course modified by the presence of 
significant local velocity dispersions, which yield profiles that 
are in some way the convolution of the above prototypes with the 
local distribution (often assumed to be described by a Gaussian 
function). A discussion on the line formation when the local 
velocity dispersion is not negligible requires detailed compu- 
tations. Therefore, we will assume in this Section that the tur- 
bulence velocity is much smaller than the macroscopic velocity 
field (which is very often the case in CSEs around AGB stars, 
with expansion velocities around 10 km s~' and local disper- 
sions ~ 1 km s~'). 

In the very inner circumstellar shells we do not expect such 
constant outwards velocities. When significant acceleration is 
still present, simple sharply-peaked profiles are expected; see 
general discussion and comparison with standard profiles in 
iBuiarrabal etai1(ll989l) . 

Simple considerations can yield analytical insight into the 
expected profiles from these regions. We recall that, in our 
sources, we can assume spherical symmetry and isotropic ex- 
pansion, thermalized populations, T k » B, unresolved sources, 
and a constant abundance. 

In the optically thin case, the velocity-integrated emission in 
a certain CO line of a certain mass of gas M at a temperature 
T k is just proportional to M/T k (from eqs. 3 to 12). Therefore, 
the brightness temperature emitted by a geometrically thin shell 
(radius: r) integrated over observable velocities, V b s (projected 
on a line of sight), would be 



T mb (r, Vobs) dV obs oc 4nr n(r)/T k 



(13) 
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Some obvious geometrical considerations show that the bright- 
ness, T mD (r, V bs), from the thin shell is independent of V b s , pro- 
vided that V bs differs from the systemic velocity, V sys , by less 
than plus/minus the expansion velocity of the shell, V(r). Since 
each shell emits at velocities between V sys - V and V sys + V, we 
can write: T mb (r, V obs ) oc 4nr 2 n(r)/T k /V(r). 

If we assume a constant mass-loss rate, M, n(r) oc 
M/r 2 /V(r), and 



T mh (r, V obs ) oc M/T k /V 2 (r) 



(14) 



The total emission of the optically thin envelope, between its 
inner and outer radii, Ri and R a , is therefore: 



7mb(Vobs) oc M l/T k (r)/V 2 (r) dr . 



(15) 



Where the integral extends only to shells such that V sys - V(r) 
< V bs - V sys + V(r), the only ones that can emit at V ODS . When 
the expansion velocity is constant, the resulting intensity does 
not depend on V bs within V sys + V, being zero otherwise; the 
profiles are therefore flat-topped, one of the cases discussed by 
lOlofsson et all (11982b . 

If we assume, for instance, a constant increase of V with r, 
V(r) = V(R )r/R , and constant T k , the variation of T mD with 
V bs from eq. 15 becomes: 

7mb(Vobs) oc „ r V(R °l : - 1 (16) 



V s , 



|V obs - 

' sysl 

Except for truncation for central velocities such that | V bs - V sys | 
< V(Ri), for which the intensity remains constant. The line shape 
then shows a central maximum, which can be quite sharp be- 
cause usually R[ is much smaller than R Q . 

If T k varies like 1/r, the solution is also sharply peaked: 

V(R a ) 



T mD (V ob s) oc In 



|V obs - V s , 



(17) 



for I V bs - V sys | > V(R{), and constant for smaller values of | V b s - 

v sys |. 

Finally, we note that the conversion of the proportionality 
relations given here into true equations (i.e. the substitution of oc 
for =) is straightforward using the complete expression for t as 
given in eq. 3. 

3.3. Line shapes: optically thick emission 

As in Sect. 3.2, our model CSE presents spherical symmetry and 
isotropic expansion. Let us assume now that the emission from a 
geometrically thin shell within a velocity range, AV b s , centered 
around the systemic velocity V sys , is optically thick. In general, 
AVobs/2 can be assumed to be equal to the maximum expansion 
velocity in the CSE. The corresponding main-beam temperature 
is then just proportional to the kinetic temperature, T k , and to 
the solid angle from which emission at this velocity is emitted 
(Sect. 3.1). Let's also assume that the expansion velocity V(r) 
is constant or monotonically increasing with the distance to the 
center, in order to simplify the formulae, and again a small local 
velocity dispersion. Also for the sake of simplicity in formulae, 
we will assume that V sys = 0; in other words, our observed ve- 
locity Vobs is considered in this subsection to be measured with 
respect to the characteristic velocity of the source. 

The total emission in main-beam temperature units of the 
circumstellar envelope (with inner and outer radii Ri and R ), 
within + AV bs/2 (i.e. within V sys + AV b s /2), is then 



r m b(Vobs) = F d 2n pT k {r) dp 



r 



(18) 



Where R v is the minimum radius emitting at V b s , taking into 
account projection on the line of sight [i.e. R v is the minimum 
r within [Ri,R Q ] such that | Vobs I ^ V(r)]. In this equation, F d 
is the inverse solid angle of the telescope beam in length units 
for a given distance to the source, D. For a Gaussian beam, 



F d (cm) 



— , with Wb(cm) = W b (arcsec)D(pc) 1.5 10 1 



and Wb(arcsec) being the half-maximum full width of the beam; 
Wb(cm) is always assumed here to be much larger than the outer 
radius R (cm), as expected in our case. Finally, p is the varying 
impact parameter, i.e. the radius in the plane of the sky of the 
differential ring that is emitting at V b s and is placed at given 
distance to the center r. Due to the projection, p satisfies: 



(19) 



Therefore, the main-beam brightness temperature in the op- 
tically thick case is given by: 



r m b(V bs) = F d 2n 



1-^L/V 2 (r) x 



V(r) - rV'(r) 



V 2 (r) 



V V2 W- y obs + 



V(r)V'(r) 



V(r) 



^VHr) - 



^obs 



dr. (20) 



These equations systematically lead to profiles with a central 
peak, but less sharp than in the optically thin case. For instance, 
if the expansion velocity i s constant, V(r) i - V e x P , we get the well 
known parabolic profiles (lOlofsson et al.lll982h : 



r m b(V b S ) = F d 27T 



1 



V 2 , v 

y obs 

V 2 , 

y exp / 



f 



rT k (r) dr 



(21) 



In the case of a constant velocity gradient, V(r) = V(R )r/R , 
which probably applies in our sources, the solution is also sim- 
ple: 



r m b(V bs) = F d 2n rT k (r) dr 



(22) 



For the assumed velocity field, R v - |Vobsl^o/V(R ); except 
when Ri > |V bsl^o/V(7? ), then R v = R[. This property leads 
systematically to a truncation of the central parts of the pro- 
file, |V bsl < V(R{), for which r m b(V b s ) is constant and equal 
to 7'mb(Vobs=V(/?i)). We note that, however, R[ is usually much 
smaller than R in circumstellar envelopes around evolved stars, 
and then the effects of truncation are negligible. 

For constant T k , parabolic profiles are again predicted: 
r m b(Vob s ) = F d R 2 nT k [1 - V 2 bs /V 2 (fl )], assuming very small 
Ri (otherwise the profiles are parabolic except for truncation in 
the central velocities). 

It is, however, more realistic to assume that the tempera- 
ture decreases outwards, which leads to more sharply peaked 
shapes. If, for instance, T k = T k (R )R /r [and, we recall, V(r) 
= V(R )r/R ], we get from eq. 22 sharp profiles that are exactly 
triangular if R[ <zzR : 



7mb(V obs ) = F d 2nT k (R )R 2 



1 



I Vol 



V(R ) 



(23) 



When the temperature decreases linearly from Ri to R and 
again for R, <k R q , eq. 22 gives: 

7mb(V bs) = F d nR 2 J k {Ri)[\ - V 2 ob JV 2 {R )] - 



-F d n[T k (Ri) - T k (R )][l ~ W ob s\ 3 /V 3 (R )] . 



(24) 
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In this case, the profile is more sharply peaked than the parabolic 
profile obtained for constant 7k, but without the discontinu- 
ity in its slope at V b s = characteristic of triangular profiles. 
Truncation also appears for |V b s | < V(^i) and would be notice- 
able in the unexpected case that R\ is not much smaller than R Q . 

These equations are particularly interesting for us, because 
probably the kinetics of the circumstellar shells in our sources 
shows a significant velocity gradient and our 12 CO lines are op- 
tically thick. We therefore expect more or less triangular pro- 
files with a central peak. As we will see below, other phenomena 
should also be taken into account to explain the asymmetry in 
the observed profiles. 

3.3.1. Selfabsorption in optically thick envelopes 

In discussion in Sects. 3.2 and 3.3, we have assumed that the lo- 
cal velocity dispersion, presumably due to turbulent movements, 
is much smaller than the macroscopic velocities. A variety of 
line profiles were found, but, in particular, in all cases the pro- 
files were symmetric around the systemic velocity V sys . If the lo- 
cal velocity dispersion is not negligible (and the temperature is 
not constant), selfabsorption must be taken into account. In the 
most usual case, in which the temperature decreases outwards, 
selfabsorption is important for the most negative velocities (gas 
approaching the observer), since then outer cool regions absorb 
the emission of hotter, inner regions. This is not the case of the 
positive velocities, in which the hotter regions would be closer 
to the observer, its emission dominating the observed line. The 
result is an asymmetric profile in which the maximum is shifted 
to positive velocities, while keeping the total emission range. 

Unfortunately, these effects are too complex to be discussed 
by means of analytical, simple formulae, as those given before, 
and numerical calculations must be performed. In Sect. 3.4 we 
will present a numerical model that takes into account selfab- 
sorption, within an accurate description of radiative transfer in 
expanding envelopes. We will also see in Sect. 4 how these phe- 
nomena are useful to better understand our observations of SSs, 
particularly the asymmetric profile of the 7=2-1 line observed 
in R Aqr. 

3.4. Numerical model of line formation 

We have used a numerical model to simu late the line fo r matio n 
in our case similar to that described by iTevssier et alj (|2006), 
adapted to the high densities and temperatures expected in the in- 
ner circumstellar shells. See more details on the numerical code 
in that paper. 

In the calculation of the level population, we have taken into 
account collisional and radiative excitation, considering a high 
number of rotational levels in two vibrational states. Our calcu- 
lations confirm that, as shown in Sect. 3.1, the population of all 
relevant levels is very accurately thermalized due to the expected 
high densities. Therefore, our results practically do not depend 
on the details of the treatment of the excitation. The radiative 
transfer equation is accurately solved taking into account a local 
turbulence velocity and the macroscopic velocity field. We so 
calculate the brightness in the direction to the observer along a 
number of lines of sight and for a number of projected velocities. 
These results are convolved with the beam shape assumed for 
the observed lines, obtaining a Rayleigh-Jeans-equivalent main- 
beam temperature for a number of LSR velocities, directly com- 
parable to the observed profiles. Thanks to the small extent of 
the source, this convolution practically does not depend on the 



exact beam shape; the beam width is therefore the only relevant 
parameter in the calculation. 

In this model, we are assuming that the emitting region is 
spherical and expands isotropically. We also assume that the 
physical and chemical conditions are in our case similar to 
those typical of inner shells around standard AGB stars. This 
may be not true in symbiotic systems, whose inner circumstel- 
lar layers could present an axis of symmetry, anomalous den- 
sity and velocity distributions, and other peculiar properties. We 
think that our lack of knowledge on these probably complex 
regions (from molecular line observations or from other stud- 
ies) prevents a more detailed modelling. We note that observa- 
tions of SiO maser emission in R Aqr ( e.g. ICotton et al.l f2004, 
iKamohara et"ai1l2010l IPardo et al.ll2004[) show properties com- 
pletely similar to those found in standard AGB stars. VLBI maps 
show, in particular, that the SiO maser in general occupy a ring- 
like region (with radius ~ 5 10 13 cm) that is comparable to those 
typical in AGB stars. Kamohara et al. found that the spatial dis- 
tribution tends, nevertheless, to show a certain axial symmetry. 
The departures from spherical symmetry do not always show the 
same pattern and are in many epochs quite small, but the axis 
orientation has remained stable at a position angle of about -10° 
during more than 10 years. The SiO data therefore indicate that 
the inner shells around R Aqr are more or less spherical and, 
because SiO masers are excited normally, must show physical 
conditions similar to those usually present in AGB stars; at least 
up to a distance of about 5 10 13 cm, i.e. in regions not much 
smaller that the total region assumed to be responsible for the 
CO lines in our model. 

4. Results 

We have detected 12 CO 7=2-1 and 7=1-0 emission from the 
symbiotic stellar systems (SSs) R Aqr and CH Cyg (the detec- 
tion of 7=1-0 in R Aqr being tentative). We also present lim- 
its of the emission of I2 CO 7=2-1 and 7=1-0 in another SS, 
HMSge, and of the emission of 13 CO 7=2-1 and 7=1-0 and 
SiO 7=5-4 in R Aqr. A summary of our observational results 
can be seen in Table 1 and the detected lines are shown in Figs. 
1, 2. The interpretation of these data can be done by means of 
the relatively simple formulation presented in Sects. 3.1 - 3.3; 
but, as we will see, some observational features require numer- 
ical calculations, using the code presented in Sect. 3.4. We will 
see that the observations are reasonably well explained by those 
predictions, and other cases that do not seem compatible with 
our data, like the narrow two-horn profiles from rotating disks, 
will not be considered. 

4.1. CO lines in R Aqr 

Let us apply the simple formulae in Sect. 3 to our observations of 
CO in R Aqr. We will first assume that the radius of the emitting 
region is smaller than the component separation, i.e. that L ~ 
0709 (3 x 10 14 cm at a distance of 214 pc). Then, the upper 
limit to the antenna temperature is given by the optically thick 
case emission: 7 m b( 12 C07=2-l) ~ 7\ x 6 10~ 5 K. The observed 
intensity in R Aqr implies a typical kinetic temperature of about 
500 K; slightly higher temperatures will be derived if the line 
opacity is not extremely large, which, as we will argue below, 
is probably the case. From eq. 3, Sect. 3.1, and assuming a CO 
abundance of about 5 x 10~ 4 , typical in CSEs around O-rich 
Mira-type stars, we deduce that the assumption of optically thick 
I2 CO 7=2-1 emission implies densities > 5 10 8 ctrT 3 . As we 
have discussed in Sect. 1 , both deduced values of the temperature 
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Table 1. Summary of observational results: peak main-beam 
temperature or 3<x limit, lcr noise, and spectral resolution with 
which the noise was calculated. 



source 


line 


peak T m b or 3<x limit 


lcr noise (res.) 






mK 


mK 


R Aqr 


12 CO 2-1 


30 


2 (2.6 km s" 1 ) 




12 CO 1-0 


~6 


2 (5.2 km s" 1 ) 




13 CO 2-1 


< 18 


6 (2.7 km s" 1 ) 




13 CO 1-0 


< 18 


6 (2.6 km s- 1 ) 




SiO 5^1 


< 19 


6 (2.8 km s" 1 ) 


CH Cyg 


12 CO 2-1 


9 


2 (2.6 km s- 1 ) 




12 CO 1-0 


-4.5 


1.4 (5.4 km s" 1 ) 


HMSge 


12 CO 2-1 


< 19 


6.5 (2.6 km s- 1 ) 




12 CO 1-0 


<26 


9 (2.6 km s- 1 ) 



V(km/s) T,(K) 




- 1000 



- 500 



5 10' 3 10 14 1.5 10 14 2 10 14 r(cm) 

Fig. 3. Distributions of the expansion velocity (left scale) and 
kinetic temperature (red dash line, right scale) deduced from our 
simple modelization of the CO-rich, inner shells around R Aqr. 



and density agree with current expectations for the innermost 
shells around AGB stars. 

We have discussed in Sect. 3.1, eqs. 6, 10, 11, and 12, 
the expected line ratios in our case. In the completely opaque 
case, the intensity ratio of the 7=2-1 and 7=1-0 lines must be 
r mb (/=2-l)/r mb (7=l-0) ~ W 2 (7=l-0)/W 2 (7=2-l) ~ 3.5. In 
the fully optically thin case, T mb (2-l)/T mb (l-0) ~ r(2-l)/r(l- 
0) x W 2 (l-0)/W 2 (2-l) = 4 w2(l-0)/W 2 (2-l) ~ 14. The line T mb 
ratio is ~ 6 (if we believe that our CO 7=1-0 line is detected, i 6 
otherwise; in any case, the uncertainty is very large). Therefore, 
our observations suggest that the CO 7=1-0 line in this source 
is moderately optically thin, typical values of t(7=1-0) ~ 0.8 
- 1 are enough to explain the (tentatively) measured line ratio. 
Since the typical r ratio is t(7=2-1)/t(7=1-0) ~ 4, the 7=2-1 
line would be in any case opaque. These values of t(7=1-0) 
imply densities ~ 10 9 crrT 3 (from formulae in Sect. 3.1), very 
reasonable in this context. The total mass in the emitting shell is 
therefore deduced to be ~ 2.5 x 10~ 5 M Q . 

Other possibilities that could explain the observations are 
much less probable. A significantly larger extent would imply 
temperature almost decreasing with 1/L 2 , to keep the observed 
intensity of the 7=2-1 line. For instance, L ~ 10 15 cm would 
imply an average temperature of ~ 45 K within the emitting re- 
gion, surprisingly low for shells around Mira-type stars (mostly 
taken into account the presence of a hot companion). Values of 
L as small as ~ 10 14 cm would imply typical temperatures of 
about 5000 K, even higher than the surface temperature of the 
cool star. 

These intensity predictions, deduced from a simple but ro- 
bust line formation model, have been fully confirmed by com- 
putations performed with the model presented in Sect. 3.4. 
(Calculations performed with this code will be discussed in more 
detail below.) 

The shape of our 12 CO 7=2-1 profile is clearly differ- 
ent from the wide profiles observed in the extended envelopes 
around standard AGB stars (parabolic, flat-topped or showing 
two horns). The single-peak, relatively sharp profile is compati- 
ble with the expected profile for the compact, very small shells 
expected in SSs (Sects. 3.2, 3.3), in which a significant acceler- 
ation of the gas is expected to be present. 

Therefore, both the intensity of the CO 7=2-1 and 7=1-0 
lines and the shape of 7=2-1 are compatible with our expecta- 



tions for the CO emission in circumstellar regions closer than 
about 2 x 10 14 cm. 

To better understand the properties of the emitting region and 
check the above conclusions, we have tried to fit our CO ob- 
servations of R Aqr with the predictions of the numerical line 
formation model presented in Sect. 3.4. We must keep in mind 
that the source modeling and the determination of the properties 
of this shell can be only a first step in the study of these lay- 
ers, because of the lack of information still existing on them, 
including, in particular, the poor data on molecular emission. 
Moreover, in the inner circumstellar shells around the cool com- 
ponent of an interacting binary system, we can expect complex 
structures and velocity fields, perhaps difficult to compare with 
the smooth, continuous ejection of material by standard AGB 
stars, and yielding line shapes with features that cannot be ac- 
counted for by our simple shell model. The best we can do is 
to find a set of characteristic properties of the emitting region, 
namely, extent, velocity field and physical conditions, that yield 
predicted line profiles compatible with the observations. 

The envelope extent is defined by its inner and outer radii, 
R„ R . We have imposed a small value of slightly larger than 
the stellar radius, and a value of R a equal to 2 x 10 14 cm, to 
be in agreement with general considerations discussed above. 
We have assumed that the velocity, V(r), increases from a cer- 
tain point in the inner envelope (supposed to be related to dust 
formation). We assume a turbulence velocity comparable to the 
minimum expansion velocity, V(Ri). The temperature is assumed 
to vary linearly in the considered region, between T^(Ri) and 
Tk{Ro), and the density varies assuming a constant mass-loss 

rate, M, i.e.: n(r) — 4m j^. (f) ■ For both the kinetic temperature and 
the expansion velocity, only linear functions are considered, in 
view of our poor knowledge on these parameters in our case. 
The CO abundance, X(CO), is assumed to be constant and equal 
to 5 x 10 , a typical value in CSEs around O-rich AGB stars. 
Note that, since the level population is practically thermalized, 
we can vary the value of X(CO) and obtain practically the same 
results, provided that the density varies in the opposite sense and 
the product nxX remains constant. 

We have found a set of values for the above parameters that 
can explain the observations. The temperature and velocity laws 
are given in Fig. 3. We note that in the outer shells around stan- 
dard AGB stars the velocity gradient becomes very small and the 
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expansion velocity increases slowly up to a certain asymptotic 
value, but this regime apparently is not reached in the compact 
shells around R Aqr we are describing here. The best-fit value 
of the mass-loss rate is 9 x 1CT 6 M yr _1 . In Fig. 4 we show the 
predicted profile, superimposed to the observed one. 

Because of the poor observational data on I2 CO 7=1-0, we 
have not tried to fit the observed profiles, but we have checked 
that the predicted intensity, T m \, ~ 5 mK, is compatible with the 
observation. 

Note that the asymmetric line shape, with a slightly red- 
shifted maximum, is well reproduced by our model, due to ab- 
sorption of the emission of inner shells by the cooler outer 
regions, only noticeable at relatively negative velocities (Sect. 
3.3.1). The expansion velocity (a parameter almost directly 
given by the observed profile width) is also quite high, compared 
to other AGB shells, compatible with the idea that mass ejection 
is particularly efficient in SSs (Sect. 1). 

Given the lack of information about the emitting region and 
the relatively simple model we are using, it is obvious that the 
values deduced here for the different parameters are indicative 
and relatively uncertain. Nevertheless, we think that the fact that 
we are able to reasonably fit the observations shows that our 
basic assumptions about the emitting region (size, temperature 
and velocity ranges, densities, ...) are essentially correct. This 
conclusion is supported by observational results of SiO maser 
emission in this object, see Sect. 3.4, which suggest that in- 
ner shells around R Aqr (not much smaller than our CO emit- 
ting region) are more or less spherical and show physical condi- 
tions similar to those typical in AGB stars. On the other hand, 
moderately oblate density distributions have been proposed for 
the inner circumstellar regions in S Ss fcawrvszc zak et al.l l2002, 
iPodsiadlow ski & Mohamed 2007), as due to gravitational fo- 
cusing. From the existing data, we cannot discriminate between 
these oblate distributions and more spherical ones. Since the bi- 
nary orbit in R Aqr is nearly edge-on, the assumption of such 
an oblate structure would limit our requirements of relatively 
high column densities to its equator. Therefore, the values of 
the mass-loss rates for R Aqr would be, in this case, somewhat 
smaller that those given above; we deduce smaller values of M 
by a factor of ~ 2. A similar smaller value of the mass-loss rate 
can be obtained, in general, if we allow the radius of the spheri- 
cal shell to increase moderately, for instance up to the separation 
between the stars, 2.3 - 2.5 10 14 cm. The data could be also 
reproduced assuming the source to be somewhat smaller than 
in our standard model, the gas temperature and the mass-loss 
rate being then deduced to be larger. (Important variations of the 
source size would lead to unexpected properties of the emitting 
gas, as mentioned before.) Strong variations in the total mass 
and mass-loss rate in the case of a significantly clumpy medium 
are also improbable, provided that the total size does not vary a 
lot, because the typical opacity and column density required to 
explain the observations would not change. 

Published estimates of the mass-loss rates from the Mira 
component of R Aqr (and in general for symbiotic Miras) are 
based either on radio continuum observations or FIR IRAS data 
(e.g. iKenvon et al.ll 19881 [Seaquist & Tavlorlll990h . The former 
method probes the ionized portion of the wind, whereas the 
later measures the amou nt of dust in the system. In particular, 
Seaq uist & Tavlorl d 19901) estimated a mass-loss rate from the 
Mira of 3.7 x 10~ 8 M yr~' from 6 cm radio flux, which would 
become 3 10 -8 M yr _1 for the distance value adopted here and 
an expansion velocity y exp ~ 15 km s~'(note that there is a typ- 
ing error in Table IV of Kenyon et al. 1988, they should give 




-50 

Fig. 4. 12 CO J -2-1 profile predicted by our best-fit model of the 
inner shells around R Aqr, red dashed line, superimposed to our 
observed profile. 



the same value for the gas mass-loss rate as Seaquist & Taylor 
because both papers used the same data). 

The IRAS data f o r sym bi otic stars were analyz ed by 
l Anandarao & Pottaschl (Il986l) . lAnandarao et alJ ([1988), and 
IKenvon et al.l d!988f) . The first group, based on analysis of the 
IRAS LRS data and the broad-band IRAS photometry, found 
two dust shells in R Aqr: a hotter inner shell with a temperature 
of about 800 K, a radius of about 6 AU or ~ 10 14 cm (assuming 
D ~ 214 pc), and a dust mass of ~ 2.5 x 10~ 7 M , which may 
overlap with the CO region, and an outer cold shell with about 87 
K, a radius of about 123 AU, and a mass of ~ 5.7 x 10~ 4 M Q . The 
hot inner shell alone can account for the LRS fluxes, whereas the 
second outer shell explains the far-IR fluxes, mainly in the 60 
and 100 micron IRAS bands. A similar complex dust structure 
was found in all of seven studied symbiotic systems, and was 
interpreted as a signature of multi ple circumstellar shells d ue to 
a discontinuous mass distribution. lAnandarao et al] {l988) also 
suggested that the inner shells represent the basic envelopes of 
the cool components, whereas the outer shells are circumbinary 
components. If the inner shells indeed result from the present 
mass ejection by the Mira, this value of dust mass requires a 
mass-loss rate of about 10 M yr _1 (assuming gas/dust mass 
ratio ~ 100 and V exp ~ 15 km s _1 ), in good agreement with our 
estimate based on CO line emission. 

On the other hand, the dust shell parameters der ived from 
broad-band IRAS photometry by IKenvon et alJ d 1 988b . temper- 
ature ~ 450 K and radius of ~ 70 AU, are significantly different 
from the above values, and they indeed are incompatible with 
the LRS data. Their dust shell mass, ~ 1.6 x 10~ 7 M (recal- 
culated for our distance), is similar to the mass of the inner 
shell proposed by Anandarao et al.; but, due to the larger ex- 
tent of the shell, the mass-loss rate is ~ 7 x 10 -7 M yr~' (as- 
suming gas/dust mass ratio ~ 100, and V exp ~ 15 kms" 1 ), an 
order of magnitude low er than our CO-based estimate. Finally, 
iGromadzki et a l. (2009) derived recently a mass-loss rate M ~ 
10~ 6 M Q yr" 1 for R Aqr, from the K-[12] color versus M rela- 
tion calibrated for normal AGB stars, but its extrapolation to our 
case is uncertain. 

In general, the comparison between the mass-loss rates de- 
rived from dust and molecular line emissions in R Aqr is very 
uncertain, because the interpretation of both molecular and FIR 
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data under our extreme conditions is not standard and because 
the dust-to-gas ratio may be significantly different in shells 
around SSs and in standard AGB stars. The dust content in the 
molecule-rich shells in symbiotic systems may be lower than in 
CSEs around normal evolved stars, because we only detect CO 
line emission from very inner shells in which grains are not yet 
completely formed, as in fact suggested by the strong velocity 
gradient deduced from our modelization. On the other hand, it is 
clear than grains can survive much more easily to the radiation 
from the hot companion, therefore we could detect dust emission 
from regions in which molecules do not exist. 

The comparison of results on ionized and molecular gas also 
needs some discussion. The high-density gas is, in general, much 
harder to ionize and photodissociate by the hot companion be- 
cause of self-shielding. We have mentioned that hydrodynam- 
ical numerical simulations show that flows from the Mira in a 
SSs can be concentrated towards the orbital plane, resulting in 
a large-scale density enhance ment in the orbital plane and low- 
density polar regions (see e.g. lGawrvszczak et al.ll2.Q02l and ref- 
erences therein). One can then expect the molecular material 
(and to a lesser extent the dust) to be placed near the orbital 
plane, as well as the presence of two ionized regions in the poles 
(and probably outside the orbit). The methods based on observa- 
tions of molecular gas and dust would refer to the denser mate- 
rial, probably lying close to the orbital plane, whereas the radio 
continuum data would probe the low-mass regions. As a result, 
one should get higher mass values from molecular lines and dust 
than from radio continuum data. Our rates from CO are about 2 
orders of magnitude higher than published values based on radio 
data, in agreement with these expectations. 

4.2. CO lines in CH Cyg 

The poor 12 CO profiles observed in CH Cyg prevent any de- 
tailed fitting of the line shape. Both lines show roughly a central 
peak and could be compatible with the emission of a region with 
significant velocity gradient. We can derive some characteristics 
of the emitting region from the line intensity. First of all, we can 
see that the J=2— 1/7= 1-0 intensity ratio is relatively high, com- 
patible with the optically thick ratio, ~ 3.5 (Sect. 3.1; see also 
discussion for R Aqr in 4.1). We also note that the component 
separation in the object is ~ 9 AU, smaller than that of R Aqr. 
Therefore we can assume that both lines are optically thick and 
come from a very compact region. From our discussion in Sect. 
3.1, we can deduce that the observed line intensities are compat- 
ible with an emitting region size (typical diameter) of about 10 
AU (typical radius ~ 5 AU, somewhat smaller than the compo- 
nent separation), and typical kinetic temperature of about 800 K. 
We stress that, because of the lack of good CO profiles in this 
source, the assumption that the emission comes from a small 
shell in expansion is less well founded than for R Aqr; this un- 
certainty and the low S/N ratio obviously result in less accurate 
conclusions from the data. 

Note that the detected profiles are, within the uncertainties, 
quite wide, suggesting that high expansion velocities are present 
in the CO-rich shells, ~ 25 km s , even larger than those mea- 
sured for R Aqr. The mass-loss rates must also be quite high, 
also larger than for R Aqr, to take into account the larger opac- 
ities and velocity dispersion (which enters both in the opacity, 
eq. 4, and in the estimate of the shell lifetime); we deduce, fol- 
lowing the prescriptions in Sect. 3.1, M t 2 10 5 M Q yr _1 . We 
suggest that the high values of the mass-loss rate and expansion 
velocity are partially due to the fact that the CH Cyg system is 
tighter than that of R Aqr. 



The high mass-loss rate in CH Cyg may be surprising in par- 
ticular because it is a semiregular (SR) variable, and it is well 
known that SR variables in general show significantly lower 
mass-loss rates than Mira-type variables. However, the case of 
CH Cyg is peculiar. CH Cyg presents a complex variation pat- 
tern that probably includes two periods of about 100 and 760 
days, plus other long term variations, with a high overall vari- 
ability amplitude, about 3 mag in the visible; see the light curve 
by the AAVSO and Mikolajewski et al. (1992). Mikolajewski et 
al. indeed classify this source as SRa variable. The basic proper- 
ties and evolutionary status of the SR varia bles of type SRa and 
SRb were compreh ensively discussed by iKerschbaum & Hronl 
(11992111994 fl996). who found that SRa stars appear as inter- 
mediate objects between Miras and SRb variables in all aspects, 
including periods, amplitudes, and mass-loss rates. Gromadzki 
et al. (2007) found that most giants in symbiotic systems reveal 
more or less regular pulsations with periods in the range 50-400 
days. They also concluded that the presence of such a variability 
can account for the relatively high mass-loss rates usually found 
in symbiotic stars as compared with single field giants (Sects. 1, 
5). 

It is also remarkable that semiregular variables with ampli- 
tudes Z 2.5 mag (like W Hya, GY Aql, T Ari, etc) tend to present 
strong SiO maser em ission, comparable to that of standard Miras 
(Alcolea et al. 1990). Therefore, the density and general physi- 
cal conditions in the inner circumstellar layers should not signif- 
icantly differ from those of Mira-type variables. These and other 
SR stars (W Hya, GY Aql, RX Boo, EP Aqr, X Her, ...) show 
dense extended envelopes well detected in CO. It is also well 
known that a number of SRs (EP Aqr, RX Boo, X Her, RS Cnc, 
etc) show extend ed shells with a clear axis of symmetry (e.g. 
iNakashimal 120051) . similar to those found in Mira-type stars in 
binary systems, like o Cet and V Hya (e.g. lKahane et alJ [T996). 

On the other hand, see Sects. 4.1, 5, it is obvious that the 
mass ejection by the evolved component of a SS can be seri- 
ously affected by the stellar interaction, and we have seen that 
mass-loss rates in such systems tend to be larger than in isolated 
AGB stars (Sects. 1, 5). It is then to be expected that the interact- 
ing nature of the CH Cyg system strongly affects the structure, 
dynamics, and density of the inner circumstellar layers, helping 
to understand the measured high amount of gas in these regions. 

The presence of a relatively high amount of mass in the in- 
ner circumstellar shells of CH Cyg is confirmed by its FIR dust 
emission and by the identification of a hot dust-shell, a few times 
larger than the stellar photosphere, that sig nificantly contributes 
to the total NIR flux JPedretti et all 120091) . The values derived 
here for the mass-loss rate and typical temperature in the inner 
shells are quite similar to those found by Taranova & Shenavrin 
(2007) from analysis of the dust FIR emission of recently ejected 
material. Our mass-loss rates are also compatible with the to- 
tal dust mass derived by Kenyon et al. (1988) and Hinkle et al. 
(2009), if we assume that dust emission comes from inner shells 
not much larger than those we are detecting in CO emission. 



4.3. Other molecular lines 

The other molecular observations in R Aqr, CH Cyg and HM Sge 
did not yield detections. The upper limits obtained for 13 CO lines 
are compatible with 12 CO/ 13 CO abundance ratios larger than 10, 
as usually found in similar objects. The nondetection of 12 CO 
to a limit of r m b(^=2-l) < 0.02 K in HM Sge, a source placed 
at more than 1 kpc, is to be expected if the properties of the 
emitting region are similar to those found for R Aqr and CH Cyg. 
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The nondetection of SiO 7=5-4 thermal (v=0) emission in 
R Aqr is more significant. This line is observed to be much 
weaker than 12 CO J =2-1. Most silicon is expected to be in gas- 
phase SiO in these O-rich regions in which dust grains are not 
formed, leading to relative abundances of about 1 - 5 x 10~ 5 . 
The Einstein A coefficient of SiO 7=6-5 is 5 x 1(T 4 s" 1 , al- 
most 1000 times larger than that of 12 CO 7=2-1. Therefore, 
the SiO 7=5-4 line should be optically thick in R Aqr, if it 
comes from a region similar to that of CO lines. However, SiO 
7=5-4 is more than 2 times weaker than CO 7=2-1, so the 
size of the SiO emitting region in R Aqr should be at least 2 
times smaller than for CO (taking into account the increase of 
temperatures in inner regions). We therefore conclude that the 
SiO thermal emission in R Aqr only can come from a relatively 
small region, probably not much larger than that measured for 
SiO masers (~ 10 14 cm in diameter, Sect. 1). This conclusion 
is interesting because SiO in standard AGB stars must be very 
abundant in more extended layers, at least in regions in which 
the grains are not completely formed, from which its intense 
mm-wave emission comes. Indeed, the extent of the SiO ther- 
mal emission i n standard AGB stars has been found to be larger 
than 10 15 cm dLucas etalJll992h . The relatively small SiO emit- 
ting region in R Aqr should therefore be due to that the molecular 
emission extent in this source is limited by photodissociation due 
to UV emission from the hot companion, rather than by dynam- 
ical distortion of the emitting region (Sect. 1, 5), because that 
phenomenon affects significantly more the SiO abundance than 
that of CO. We recall that CO is the most extended molecule 
in CSEs because self-shielding yields less efficient photodisso- 
ciation for this mole c ule than for less abund ant species, see e.g. 
iMamon et all (ll988h . lwiTiacv & Miliar! (Il997h . 

5. Conclusions 

We present observations of the 12 CO and 13 CO 7=2-1 and 
7=1-0 transitions and of the SiO 7=5-4 one in three symbi- 
otic stellar systems (SSs), R Aqr, CH Cyg, and HM Sge. 12 CO 
7=2-1 and 7=1-0 emissions were detected in R Aqr and CH 
Cyg. An accurate line profile of 12 CO 7=2-1 in R Aqr was ob- 
tained. 

The observed lines are very weak. If we compare them with 
those usually obse rved in standard evolv ed stars with similar 
properties (see e.g. Bujarraba l et alJ[T992h . the CO lines in SSs 
are ~ 100 times weaker, even if we also consider objects, like 
young planetary nebulae, in which molecule photodissociation 
seems to be already efficient. 

The weak intensities of the observed lines suggests that CO 
emission comes from a very compact region, much smaller than 
the region responsible for low-7 CO emission in circumstel- 
lar envelopes (CSEs) around standard AGB stars (typical radius 
larger than 10 16 cm). This result is compatible with observations 
(mostly nondetections) in SSs of other molecular lines, like the 
SiO, H2O and OH masers, see Sect. 1. We have shown, from 
very general considerations, Sects. 3, 4, that the detected CO in- 
tensities are coincident with those expected for emission coming 
from a circumstellar region comparable in radius to (or slightly 
smaller than) the distance between the stellar components of the 
system. We assumed that these molecule-rich layers in SSs have 
physical conditions similar to those of the inner layers around 
standard AGB stars. Particularly satisfactory is the fitting of the 
lines observed in the most intense source, R Aqr, in which the 
stars are separated by typically 2.3 x 10 14 cm, see Sect. 1. 

We have also performed model calculations of the CO emis- 
sion of these shells, assuming emitting region radii of that or- 



der and the physical conditions and dynamics expected in gas 
placed at similar distances around standard AGB stars (Sects. 3, 
4). Our calculations reproduce the observed intensities, within 
the observational uncertainties. We can also reproduce the 12 CO 
7=2-1 line profile observed in R Aqr, see Sect. 4.1 and Fig. 
4. Our model explains in particular the asymmetry observed in 
the profile, as a result of selfabsorption in an optically thick line 
and in the presence of significant gas acceleration and system- 
atic temperature decrease with radius. All these properties are 
expected to be present in the very inner layers from which we 
are assuming that the CO emission comes. 

The data obtained in CH Cyg are worse and their analysis 
is more uncertain. To explain the observed intensities, we de- 
duce a typical radius of about 10 14 cm, again slightly smaller 
than the distance between the stars (Sects. 1, 4.2). In this region, 
the physical conditions seem also comparable to those typical 
of AGB stars, although the expansion velocity, ~ 25 km s , is 
slightly larger than those usual in CSEs around standard red gi- 
ants. 

We so conclude that the nebulae around the SSs R Aqr and 
CH Cyg include a very compact molecule-rich region, with radii 
~ 10 14 - 2x 10 14 cm, from which the molecular emission de- 
tected in these objects comes. This region is in extent compa- 
rable to or slightly smaller than the region within both compo- 
nents of the binary system. We deduce in these regions outwards 
expansion velocities of about 5-25 km s , with a significant 
acceleration of the gas, temperatures decreasing with radius be- 
tween about 1000 and 500 K, and densities ~ 10 9 - 3 x 10 8 cm 4 . 
The general physical conditions and dynamics in this inner shell 
would be then similar to those typical of standard AGB stars. 
This result is supported by the 'normal' observational properties 
of SiO maser emission from R Aqr (Sects. 3.4, 4), which show 
that the CSE around this star is relatively normal (i.e. that its 
main properties are comparable to those of standard AGB stars) 
up to distances ~ 5 10 13 cm, at least. 

We find relatively high values of the gas mass and mass- 
loss rate. In R Aqr, we estimate that the total mass of the CO 
emitting region is ~ 2 - 3 x 10~ 5 M , and that the value of 
M ranges between 5 10~ 6 and 10~ 5 M yr _1 , in agreement with 
results obtained from dust emission. These relatively high val- 
ues therefore confirm the general trend of SSs to present sig- 
nificantly high er mass-loss rates than standard, isolat ed evolved 
stars (Mikolaiewska 1999, Mikolaiewska et al . 2002). We must 
take into account, however, that the comparison between mass- 
loss rates in the very extended and isotropic winds from standard 
AGB stars and the values of M deduced for our sources is diffi- 
cult. Here, we are really detecting gas which is confined within 
the orbit of the companion, and its future dynamics is difficult 
to know. It is probable that most of this material will never leave 
the system (in this case we would not be dealing with a true mass 
ejection) or will leave it by means of outbursts or bipolar jets. It 
is also possible that the process of mass ejection in these stars 
is significantly affected by the companion and cannot be easily 
compared with that in isolated AGB stars. Hydrodynamical sim- 
ulatio ns dGawrvszczak et all 120021 iPod siadlowski & Mohamed 
2007) show indeed that, in systems with orbital parameters sim- 
ilar to those of CH Cyg and R Aqr, the radius of the dust shell is 
comparable with the Roche lobe radius and a wind Roche-lobe 
overflow will occur. 

In any case, we think we can safely conclude that the mass 
of the gas in the close surroundings of R Aqr and CH Cyg is 
relatively high, and that its presence is probably related to the 
gravitational effects of the nearby companion, even if these pro- 
cesses are not well understood. 
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We note that the few symbiotic systems showing molecular 
emission (as R Aqr and CH Cyg) are objects showing pa rticu- 
larly large distances between the stars (Sc hwarz et al.lll995l) . It 
is then probable that in other SSs the molecule-rich region is an 
extremely small shell tightly surrounding the Mira-type compo- 
nent, which would explain its very weak molecular emission. 

From the properties we derive for the CO-rich region and 
the nondetection of SiO thermal lines (Sect. 4.3), whose emit- 
ting region must be significantly smaller than for CO, we sug- 
gest that the small extent of the molecule-rich gas in these 
sources is mainly due to molecule photodissociation by the hot 
dwarf star, which must be significantly more efficient for low- 
abundance molecules. Ionized gas has been found to be very 
ext ended in SSs, particularly in the well s tudied R Aqr (see 
e.g. lCorradi et alj2003llHollis et al.lll999albl) . where bipolar and 
disk-like ionized-gas structures are known to extend over more 
than one arcminute (several 10 17 cm), apparently without any 
molecule-rich counterpart. It is obvious that the relatively stan- 
dard inner shells around the red giant show very different prop- 
erties than the spectacular outer parts of the nebula. Therefore, 
disruption of the CSE due to the strong stellar interaction in SSs 
must also play an important role in the confinement of gas (in 
particular, molecule -rich gas) in a very compact and dense shell. 
Seaq uist et al.l d 19951) proposed that the powerful wind from the 
companion may sweep much of the material ejected by the Mira 
star. This effect would explain the strong discontinuity found be- 
tween the general properties of the region s inside and outsid e 
the orbit of the companion and, following Seaquist et al] d 19951) . 
could increase the efficiency of photodissociation in the rest of 
the nebula, helping to explain the lack of molecular emission in 
SSs. 
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